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INTRODUCTION 


This  report  summarizes  the  activities  performed  by  Molecular  Physics  Laboratory  personnel 
under  SRI  Project  P 1 7002  (Subcontract  No.  0000 1131,  The  Trustees  of  Princeton  University).  The 
Principal  Investigator  for  this  project  was  Dr.  Jochen  Marschall  and  the  Co-Investigator  was  Dr. 

Dusan  Pejakovic.  The  period  of  performance  was  1 1/4/2005  to  3/31/2009. 

This  SRI  subcontract  supports  a  prime  contract  awarded  to  Princeton  University  by  the  Air 
Force  Office  of  Scientific  Research  entitled,  “Characterization  of  Near  Wall  Surface  Chemistry  and 
Fluid  Interaction  in  Hypersonic  Boundary  Layers”,  with  Prof.  Maria  P.  Martin  as  the  Principal 
Investigator.  The  statement  of  work  for  the  SRI  subcontract  is: 

SRI  will  use  a  h  igh-power  microwave  discharge  flow  tube  coupled  with  mass  spectroscopic 
and  laser-based  diagnostics  to  investigate  heterogeneous  processes  on  heated  test  specimens  as  a 
function  of  sample  temperature,  gas  pressure,  and  gas  composition.  The  microwave  discharge  facility 
will  enable  testing  in  clean,  controlled  gas  environments  with  high  levels  of  atomic  oxygen  and  atomic 
nitrogen.  Silicon  carbide,  carbon  and  other  specimens  will  be  tested,  in  quantities  and  with 
specifications  to  be  determined  by  mutual  discussion. 

SRI  will  use  its  expertise  in  optical,  laser-based,  and  mass  spectroscopic  techniques  to  quantify 
the  chemical  composition  of  the  test  environment,  and  to  monitor  changes  in  reactant  and  product 
concentrations  as  a  function  of  various  experimental  parameters.  Changes  in  sample  mass  and 
morphology  will  be  documented  and  correlated  with  gas-phase  diagnostics.  The  experimental  results 
will  be  used  to  construct  finite-rate  surface  chemistry  models  for  use  in  computational  fluid  dynamics 
simulations. 

During  the  course  of  this  study,  the  decision  was  made  to  focus  on  catalytic  processes  on  silica 
surfaces  in  partially  dissociated  oxygen-nitrogen  gas  mixtures.  This  choice  was  motivated  by  several 
factors,  including  i)  many  reusable  thermal  protection  systems  either  have  silica-based  protective 
coatings  or  form  silica-rich  oxides  scales  upon  oxidation,  ii)  limited  experimental  information  is 
currently  available  in  the  literature  on  the  products  and  production  channels  of  O-  and  N-atom  surface 
recombination,  and  iii)  catalytic  processes  are  the  logical  starting  points  for  developing  and  testing 
finite-rate  surface  reaction  models  before  extending  such  models  to  surface  consuming  processes. 

The  following  section  presents  some  background  on  the  problem  of  O  and  N  surface 
recombination  and  its  importance  to  aerothermal  heating. 

The  next  section  describes  our  initial  set  of  experiments  which  were  aimed  at  confirming  and 
extending  previous  room-temperature  measurements  suggesting  the  importance  of  nitric  oxide  (NO)  as 
surface  recombination  product  in  partially  dissociated  oxygen-nitrogen  mixtures.  These  experiments 
used  a  diffusion  tube  side-arm  reactor  in  conjunction  with  two-photon  laser-induced  fluorescence 
(LIF)  detection  of  atomic  oxygen  and  atomic  nitrogen.  With  modeling  support  by  our  Princeton 
collaborators,  Lian  Duan  and  Pino  Martin,  our  experimental  results  demonstrate  that  NO  formation  is 
as  significant  a  pathway  as  O2  and  N2  formation.  The  results  of  these  studies  were  presented  as  a 
conference  paper1  at  the  46th  ALAA  Aerospace  Sciences  Meeting  and  Exhibit,  January  7-10,  2008  in 
Reno,  NV  and  later  published  in  the  AIAA  Journal  of  Thermophysics  and  Heat  Transfer.2 

The  last  section  of  the  report  describes  our  follow-on  experiments  aimed  at  detecting  the 
surface  production  of  NO  directly  at  elevated  temperatures,  using  single-photon  LIF  excitation  scans 
as  an  NO  diagnostic,  together  with  isotopically-enriched  l?Nl60  and  r5NftO  titration  gases  to  isolate 
surface-produced  NO  from  all  other  possible  sources.  These  experiments  independently  confirm  the 


surface  production  of  nitric  oxide  and  provide  data  that  can  be  compared  with  the  predictions  of  finite- 
rate  surface  reaction  models. 


BACKGROUND 

The  heterogeneous  recombination  of  dissociated  species  on  the  surface  of  thermal  protection 
system  (TPS)  materials  is  known  to  contribute  to  the  aerothermal  heating  of  hypersonic  reentry 
vehicles.  In  the  Earth’s  atmosphere,  oxygen  and  nitrogen  atoms  are  generated  in  the  shock  layer  ahead 
of  leading  edges  during  hypersonic  flight.  Modeling  O-  and  N-atom  surface-catalyzed  recombination 
is  an  important  part  of  computational  fluid  dynamics  (CFD)  simulations  that  predict  aerothermal 
heating  as  a  function  of  flight  trajectory. 

Traditionally,  CFD  aerothermal  heating  computations  have  treated  oxygen  and  nitrogen 
recombination  as  two  uncoupled  parallel  processes,  O  +  O  — »  O2  and  N  +  N  — >  N2  neglecting  the 
possibility  of  NO  formation  via  O  +  N  surface  reactions.  This  choice  is  partly  a  matter  of 
computational  convenience  and  partly  due  to  the  lack  of  experimental  information  on  the  importance 
of  the  NO  formation  pathway.  More  recent  modeling  efforts  have  introduced  finite  rate  kinetic  models 
to  better  capture  the  physics  of  surface  catalytic  reactions.3"9  In  at  least  one  case,  the  inclusion  of  NO 
surface  formation  in  a  finite-rate  surface  chemistry  model  seems  to  improve  agreement  between 
aerothermal  heating  computations  and  measured  flight  data.  Strong  support  for  NO  surface  formation 
is  also  provided  by  the  in  situ  laser-induced  fluorescence  (LIF)  detection  of  different  NO 
concentrations  in  front  of  different  materials  tested  in  a  plasma  wind  tunnel.1" 

Copeland  et  al.  reported  an  experimental  approach  for  inferring  the  formation  of  NO  in 
N2/N/O  gas  mixtures  using  two-photon  LIF  detection  of  N  atoms  in  a  room  temperature  diffusion-tube 
side-arm  reactor.  In  the  side-arm  reactor  technique,  atomic  species  are  generated  in  a  main  flow  tube 
and  diffuse  into  the  opening  of  a  dead-end  side-arm  tube  where  they  are  lost  by  gas-phase  and  surface- 
mediated  recombination  reactions,  establishing  a  steady-state  decaying  atom  concentration  profile 
along  the  length  of  the  side-arm  tube.  At  sufficiently  low  pressures,  heterogeneous  atom  loss 
processes  dominate  over  (most)  gas-phase  reactions,  and  this  atom  concentration  profile  can  be  related 
to  the  catalytic  activity  of  the  tube  surface.12  |Q 

Copeland  et  al.1 1  produced  a  range  N2/N/O  gas  mixtures  in  the  main  flow  tube  by  titrating  a 
partially  dissociated  nitrogen  flow  with  controlled  amounts  of  NO  to  drive  the  gas-phase  reaction 
N  +  NO  — >  N2  +  O.  This  bimolecular  reaction  produces  one  O  atom  for  every  NO  consumed  and  is 
very  rapid  with  a  room-temperature  rate  coefficient  of  3  x  10'  cm  molecule"  s"  .  They  then 
demonstrated  that  as  O  atoms  were  added  to  the  flowing  N2/N  mixture,  the  loss  of  N  atoms  from  the 
gas-phase  within  the  side-arm  tube  increased.  They  interpreted  this  accelerated  N-atom  loss  as 
resulting  from  the  same  gas  phase  N  +  NO  — >  N2  +  O  reaction,  where  the  NO  source  was  now  N  +  O 
surface  recombination  on  the  walls  of  the  side-arm  tube. 

Under  the  low  pressures  of  their  experiment  (13  to  65  Pa),  the  direct  production  of  NO  in  the 
gas  phase  by  the  three-body  N  +  O  +  M  reaction  was  insignificant.21  Because  the  N  +  NO  — »  N2  +  O 
reaction  is  so  fast  and  atom  concentrations  in  the  side-arm  are  low,  the  steady-state  gas-phase  NO 
concentrations  generated  by  surface  recombination  are  too  small  for  direct  LIF  detection,  and  NO 
surface  formation  must  be  inferred  indirectly  from  observed  changes  in  atom  concentrations. 

Copeland  et  al.1 1  suggested  a  further  check  of  their  interpretation:  if  surface  NO  formation  was 
indeed  occurring,  O-atom  loss  from  the  gas-phase  should  decrease  with  the  addition  of  N  atoms,  since 


2 


each  reaction  event  N  +  NO  regenerates  a  gas-phase  O-atom.  However,  they  did  not  perform  O-atom 
LIF  measurements  to  confirm  this. 


SURFACE  PRODUCTION  OF  NITRIC  OXIDE  AT  ROOM  TEMPERATURE:  O-  AND  N- 
ATOM  DETECTION  EXPERIMENTS 

We  have  repeated  and  extended  the  measurements  of  Copeland  et  al.  by  adding,  i)  a  more 
powerful  discharge  source  for  N-atom  production,  ii)  the  capability  for  sequential  N-atom  and  O-atom 
LIF  diagnostics,  and  iii)  simultaneous  concentration  measurements  at  four  different  side-arm 
locations. 
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Figure  1.  Experimental  setup;  PMT,  SHG,  and  THG  indicate  photomultiplier  tube,  second  harmonic 
generation  and  third  harmonic  generation,  respectively. 
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Experimental  Setup 

The  layout  of  our  experiment  is  shown  in  Fig.l .  A  partially-dissociated  nitrogen  flow  is 
produced  by  metering  nitrogen  through  a  6  kW  microwave  discharge  using  a  500  standard  cubic 
centimeter  per  minute  (seem)  Tylan  mass  flow  controller.  The  main  tube  is  made  of  glass  and  the 
side-arm  tube  is  made  of  GE214  quartz.  With  valve  1  open  and  valve  2  closed,  the  partially- 
dissociated  gas  flows  down  the  main  arm,  past  the  opening  of  the  dead-end  side-arm  tube,  and  into  a 
roots  blower  vacuum  system.  In  this  configuration  the  side-arm  operates  as  a  diffusion  tube.  The 
side-arm  can  be  reconfigured  as  a  flow  tube  by  closing  valve  1  and  opening  valve  2.  Gas  pressures  are 
measured  near  the  opening  (port  A)  and  near  the  end  (port  B)  of  the  side-arm  tube  using  a  10  Torr 
capacitance  manometer  with  four  digit  precision  and  a  reading  error  less  than  1  %. 

Species  detection  was  accomplished  using  well-known  2-photon  LIF  schemes:  excitation  at 
226  nm  followed  by  fluorescence  at  845  nm  for  O-atoms22"'  and  excitation  at  207  nm  followed  by 
fluorescence  at  745  nm  for  N-atoms.  4  For  O-atom  LIF,  a  Lambda-Physik  ENG  103  XeCl  Excimer 
laser  was  used  to  pump  a  Lambda-Physik  FL  3002  dye  laser  filled  with  Coumarin  47  laser  dye 
dissolved  in  methanol.  The  resulting  tunable  blue  light  near  452  nm  was  frequency-doubled  using  a 
type-B  beta  barium-borate  (BBO)  crystal  to  produce  226  nm  radiation.  For  N-atom  LIF,  tunable  red 
light  near  621  nm  was  produced  by  pumping  a  Quanta-Ray  MCI  dye  laser  filled  with  Rhodamine  640 
dissolved  in  methanol,  with  the  second  harmonic  output  of  a  Quanta-Ray  DCR-3  Nd:YAG  laser.  The 
red  light  was  frequency-tripled  to  produce  207  nm  radiation  using  type-2  and  type-C  BBO  crystals  for 
frequency  doubling  and  frequency  mixing,  respectively. 

The  two  ultraviolet  laser  beams  were  alternately  directed  through  a  quartz  Brewster  angle 
window  down  the  centerline  of  the  side-arm  tube  using  different  combinations  of  dichroic  mirrors, 
quartz  prisms  and  collimating  optics.  A  Pellin-Broca  prism  was  used  to  separate  visible  from 
ultraviolet  light  during  the  O-atom  detection  experiments;  dichroic  mirrors  were  sufficient  to  separate 
visible  from  UV  light  during  the  N-atom  experiments.  Quartz  flats  were  used  to  direct  small  fractions 
of  the  laser  beam  into  a  Molectron  pyroelectric  energy  meter  and  into  a  Coherent  wavemeter,  to 
monitor  laser  pulse  energy  and  wavelength,  respectively.  The  ultraviolet  excitation  energy  per  laser 
pulse  was  ~0.5  mJ  for  O-atom  detection  and  -0. 1  mJ  for  N-atom  detection. 

Fluorescence  was  collected  at  right  angles  to  the  side-arm  tube  by  four  red-sensitive 
Hamamatsu  photomultiplier  tubes  (R636)  in  gated  sockets,  fitted  with  3-nm  bandwidth  interference 
filters  centered  on  845  nm  for  O-atom  detection  and  745  nm  for  N-atom  detection.  PMT  signals  were 
amplified  by  a  factor  of  five  and  averaged  over  150  ns  wide  gates  using  Stanford  Research  Systems 
(SRS)  boxcar  integrators.  Data  acquisition  and  laser  firing  were  synchronized  at  10  Hz  using  a  SRS 
digital  delay  generator.  Automation  of  the  experiment  and  data  collection  was  performed  by  a 
personal  computer  using  LabView  software. 

The  concentration  of  O  atoms  in  the  main  gas  flow  was  controlled  by  adding  a  1  %NO  in  N2 
gas  mixture  to  the  partially-dissociate  nitrogen  flow  upstream  of  the  side-arm  opening.  At  the  N  +  NO 
titration  endpoint  all  N  atoms  are  consumed  and  replaced  in  the  flow  by  O  atoms.  The  titration 
endpoint  can  be  detected  by  the  disappearance  of  the  N-atom  LIF  signal  and  the  leveling  off  of  the  O- 
atom  LIF  signal. 

The  experiments  were  conducted  at  298  K,  with  a  44.4  seem  flow  of  N2  into  the  main  arm  and 
the  discharge  operating  at  40%  of  its  maximum  capacity  (corresponding  to  -2.4  kW.)  The  pressure  at 
port  A  before  any  addition  of  the  NO/N2  gas  mixture  at  the  titration  port  was  adjusted  to  be  60  Pa  in 
both  flow  tube  and  diffusion  tube  experiments,  by  throttling  the  flow  with  the  appropriate  valve.  The 
maximum  gas  flow  added  at  the  titration  port  was  -6  seem,  which  increased  the  pressure  at  port  A  by 
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about  1 .3  Pa  in  both  flow  and  diffusion  tube  modes.  In  flow  tube  mode,  the  pressure  difference 
between  ports  A  and  B  was  approximately  5.3  Pa  for  all  levels  of  added  titration  gas  flows,  resulting  in 
bulk  flow  velocities  ranging  from  about  3.6  to  3.85  m/s  at  the  different  PMT  locations.  Important 
experimental  dimensions  are  shown  in  Fig.  2.  The  spacing  between  PMTs  is  25.4  cm  and  the  spacing 
between  pressure  ports  A  and  B  is  120  cm.  The  inner  diameter  of  the  side-arm  tube  is  2.2  cm. 


Qualitative  Behavior 

Figures  3a-d  show  the  variation  in  relative  atom  concentration  as  a  function  of  the  flow  rate  of 
NO  added  at  the  titration  port  for  both  N  atoms  and  O  atoms  with  the  side-arm  configured  as  either  a 
flow  tube  or  a  diffusion  tube.  The  raw  LIF  data  collected  by  each  PMT  were  corrected  for  scattered 
light  contributions  by  subtracting  the  off-resonance  signals  (corrections  were  always  3%  or  less  of  the 
total  signal),  normalized  by  the  square  of  the  laser  pulse  energy,  and  averaged  over  at  least  200  laser 
shots  at  each  added  NO  flow  setting.  The  averaged  data  reflect  changes  in  atom  concentration  at  the 
particular  PMT  location  with  NO  additions  at  the  titration  port.  To  better  illustrate  common  behaviors 
in  Figs.  3a-d  the  averaged  values  were  then  put  on  a  common  scale  by  dividing  each  PMT  data  set  by 
its  maximum  value. 
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Figure  3.  Variations  in  the  scaled  atom  concentrations  measured  at  the  four  PMT  locations  with  the 
addition  of  NO  at  the  titration  port:  a)  N  atoms  in  flow  tube,  b)  N  atoms  in  diffusion  tube,  c)  O  atoms  in 
flow  tube,  and  d)  O  atoms  in  diffusion  tube.  Solid  lines  are  guides  for  the  eye;  the  offset  dashed  line  in 
Fig  3a  is  discussed  in  the  text. 


With  the  side-arm  tube  configured  as  a  flow  tube,  both  the  N-atom  and  the  O-atom 
concentrations  vary  approximately  linearly  as  NO  is  added  at  the  titration  port  (except  in  the  region 
near  zero  added  NO  in  the  N-atom  figures;  this  is  discussed  in  the  following  subsection.)  When  the  * 
reactor  is  switched  to  diffusion  tube  mode  clear  deviations  from  this  linear  dependence  are  evident. 

In  Fig.  3b  the  N-atom  concentrations  decrease  faster  with  added  NO,  producing  a  concave 
profile,  while  in  Fig.  3d  the  O-atom  profile  becomes  convex,  consistent  with  a  faster  rise  in  O-atom 
concentrations  with  added  NO.  The  result  shown  in  Fig.  3b  reproduces  the  experimental  observations 
of  Copeland  et  al.1 1  for  N  atoms,  and  the  result  shown  in  Fig.  3d  confirms  their  predictions  for  O 
atoms.  Copeland  et  al.  1  also  found  a  trend  of  increasing  deviation  from  linearity  with  increasing 
distance  from  the  side-arm  opening.  Despite  experimental  scatter,  similar  trends  are  found  in  the 
present  measurements,  and  especially  when  comparing  PMT1  (solid  circles)  and  PMT4  (open 
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triangles)  data  points  near  the  midpoint  of  the  titration  curve.  Thus,  the  current  measurements  support 
and  strengthen  the  hypothesis  of  NO  formation  by  surface  reactions. 

In  order  to  compare  experimental  data  with  the  predictions  of  a  multi-species  reaction- 
diffusion  model,  the  scaled  measurements  of  Figs.  3b  and  3d  must  corrected  for  variations  in  LIF 
detection  sensitivities  at  the  different  PMT  locations  and  the  absolute  atom  concentrations  at  the  side- 
arm  opening  must  be  determined. 


Concentration  Estimates 

Estimates  of  the  atom  concentrations  at  the  side-arm  opening  can  be  made  from  the  titration 
endpoints  under  flow  tube  conditions,  however  there  are  two  complications:  determining  the  exact 
endpoint  and  interpreting  the  initial  rise  of  the  N-atom  signal  with  the  addition  of  small  amounts  of 
NO  (Fig  3a.) 

The  titration  endpoints,  where  the  N-atom  LIF  signal  vanishes  and  the  O-atom  LIF  signal 
levels  off,  are  roughly  0.04  seem  NO  in  Fig.  3a  and  0.045  seem  NO  in  Fig.  3c.  This  difference  is 
likely  a  reflection  of  day-to-day  repeatability  (N-atom  and  O-atom  measurements  were  conducted  on 
different  days),  as  well  as  experimental  uncertainty  associated  with  detecting  a  sharp  endpoint. 

Below  the  endpoint,  the  O-atom  concentration  rises  linearly  from  zero  with  NO  addition  until 
the  endpoint  is  reached  where  all  atomic  nitrogen  has  been  consumed.  Above  the  endpoint  the  O-atom 
concentration  should  remain  constant.  If  we  choose  0.0425  seem  NO  as  the  average  endpoint,  the  O- 
atom  concentration  is  computed  from  the  measured  temperature,  pressure  and  gas  flows  as 
~2.1 1  x  10  *  mol  m '3  with  a  1-guncertainty  of  -25%.  It  is  shown  in  the  following  section,  that  O-atom 
losses  on  the  walls  between  the  titration  port  and  side-arm  entrance  are  small  and  easily  absorbed  in 
this  error  estimate. 

If  the  N  +  NO  — >  N2  +  O  reaction  is  the  dominant  reactive  process,  the  O-atom  endpoint 
concentration  should  also  equal  the  N-atom  number  density  without  added  NO.  However,  the  N-atom 
concentration  profiles  in  Figs.  3a  and  3b  show  an  initial  sharp  increase  with  the  addition  of  small 
amounts  of  NO,  in  contrast  to  the  immediate  N-atom  loss  anticipated  from  the  fast  N  +  NO  reaction. 
Similar  behaviors  have  been  documented  by  other  researchers  working  with  fast  flow  and/or  high 

25  27  ° 

power  nitrogen  discharge  systems. 

Two  main  explanations  have  been  proposed:  1)  NO  addition  causes  a  reduction  in  N-atom 
surface  losses,  or  2)  NO  addition  catalyzes  additional  dissociation  of  electronically  or  vibrationally- 
excited  nitrogen  molecules,  N2*.  It  is  possible  to  construct  kinetic  models  for  both  mechanisms  that 
can  mimic  the  experimental  observations. 25 

In  our  system,  NO  is  injected  into  the  flow  near  the  tube  centerline  and  causes  a  very  rapid  rise 
in  the  N-atom  LIF  signal  even  when  the  concentration  of  NO  is  -10  times  lower  than  the  N-atom 
concentration.  In  order  to  decrease  heterogeneous  losses  of  N  atoms  either  NO  or  the  O  atoms 
produced  by  N  +  NO  reaction  must  reach  the  surface  and  dominate  the  competition  for  surface  sites. 

It  seems  unlikely  that  NO  could  dominate  surface  coverage  under  these  conditions,  since  N  atoms  are 
both  in  excess  and  should  have  a  much  stronger  driving  force  (energy  reduction)  for  chemisorption. 
Density  functional  theory  computations  by  Arasa  et  al.  ,  suggest  O  atoms  have  a  slightly  greater 
adsorption  energy  on  Si02  than  N  atoms,  but  that  both  energies  are  high  and  initial  sticking 
coefficients  are  similar."  More  recent  atomistic  simulations  of  Cozmuta  4  predict  that  O-atom 
adsorption  energies  are  higher  than  N-atom  adsorption  energies  on  silica  surfaces  and  that  O-atom 
adsorption  dominates  over  N-atom  adsorption  in  gas  mixtures  containing  both  species. 
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The  present  experiments  were  carried  out  at  similar  total  pressures,  temperatures  and  flow  rates 
to  those  of  Copeland  et  al.1 1  who  did  not  observe  an  initial  rise  in  N-atom  concentration  with  added 
NO.  It  is  possible  that  Copeland  et  al.  did  not  collect  titration  data  with  sufficient  resolution  in  added 
NO  flow  to  observe  the  initial  N-atom  signal  rise.  However  our  microwave  discharge  source  also 
operates  at  much  higher  powers  (2.4  kW  vs.  100  W)  and  is  expected  to  generate  much  higher 
concentrations  of  excited  nitrogen  species.  Any  reduction  of  N-atom  surface  recombination  because 
of  NO  addition  should  be  similar  in  both  sets  of  experiments,  however  catalytic  dissociation  of  N2* 
upon  NO  addition  would  be  more  easily  observed  in  our  discharge  system. 

Within  the  catalytic  dissociation  model,  the  initial  rise  in  the  N-atom  concentration  is  explained 
by  a  net  reaction  NO  +  N2*— >  NO  +  N  +N,  where  the  reaction  rate  is  faster  than  the  N  +  NO  rate.  The 
subsequent  N-atom  concentration  maximum  and  a  linear  decay  with  added  NO  is  explained  by  the 
depletion  of  N2*,  the  concentration  of  which  is  presumably  lower  at  the  titration  port  than  that  of  N 
atoms.  Since  the  slope  of  the  linear  decay  then  simply  reflects  removal  of  N  atoms  by  NO,  the 
endpoint  corresponding  to  the  initial  N-atom  concentration  can  be  estimated  from  a  parallel  line  (the 
dashed  line  in  Fig.  3a)  as  approximately  0.035  seem,  or  about  12.5%  lower  than  the  observed 
endpoint.  The  initial  N-atom  concentration  corresponding  to  an  average  observed  endpoint  of  0.0425 
seem  NO  is  then  - 1 .84  x  1 0"'  mol  m'\ 


LIF  Signal  Corrections 

The  LIF  signal  magnitudes  collected  from  each  PMT  depend  on  many  factors  in  addition  to  the 
local  atom  concentration,  such  as  optical  alignment,  PMT  sensitivity,  filter  transmission,  and 
excitation  volume  and  laser  beam  geometry.  For  a  stable  species,  these  factors  can  be  determined  for 
each  PMT  by  filling  the  reactor  with  a  uniform  static  concentration  and  comparing  collected  LIF 
signals,  but  this  is  not  possible  for  reactive  species  like  O  and  N  atoms  that  are  lost  on  the  reactor 
walls.  The  closest  approximation  to  uniform  atom  concentrations  that  can  be  obtained  in  our  system  is 
under  flow  tube  conditions.  However,  the  signal  ratios  collected  in  flow-tube  mode  still  need  to  be 
adjusted  for  differences  in  atomic  concentration  at  the  different  PMT  locations  caused  by  pressure 
gradients  and  heterogeneous  atom  losses,  before  they  can  be  used  to  normalize  diffusion  tube 
measurements. 

The  Hagen-Poiseulli  relation  "  can  be  used  together  with  the  ideal  gas  law,  the  measured 
pressures  at  ports  A  and  B,  the  measured  N2  flow  rates  and  the  known  reactor  distances  to  compute  the 
total  gas  pressure  at  each  PMT  location,  assuming  that  molecular  nitrogen  dominates  the  transport 
properties  with  a  viscosity  of  17.77  x  10’*'  Pa-s  at  298  K.M  The  bulk  flow  speed,  v,  can  be  computed 
at  each  location  from  the  pressure  and  N2  flow  rate,  using  mass  conservation  and  the  ideal  gas  law. 

The  atom  loss  due  to  heterogeneous  recombination  between  PMT  locations  are  estimated  from  a 
simple  one-dimensional  first-order  model  that  equates  the  difference  in  the  atom  flow  between  z  and 
z  +  Sz  to  the  losses  incurred  on  the  wall  over  this  distance.  This  model  is  described  by  the  differential 
equation 


d(vc)  _  yv  _ 
~dz  ~~2RC  ' 


(la) 


where  R  is  the  tube  radius,  c  is  the  atom  concentration,  v  is  the  bulk  flow  speed,  v  is  the  average 
thermal  speed  of  the  atoms  ,  and  y  is  the  estimated  loss  probability  (the  fraction  of  atom- 
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surface  collisions  which  result  in  permanent  removal  of  the  atom  from  the  gas  phase).  The  value  of 
y  is  uncertain,  but  as  a  first  approximation  can  be  taken  in  the  range  -1-10  x  10  ^  for  O  and  N  atoms 
on  quartz.1 '  Defining  the  representative  atom  concentration  over  Sz  as  c  =  vc/vmt  allows  integration 
of  Eq.  (la)  to  give 


Yv{zi-zx) 

2  Rv„ 


\ 


J 


=  yexP 

r2  \ 


Mz2-zi)\ 

2  Ryave  )' 


(lb) 


where  vave is  the  average  flow  velocity  between  locations  z,  and  z2,and  v2/v,  =P]/P2  follows  from 

the  ideal  gas  law  and  mass  conservation.  Eq.  (lb)  allows  estimation  of  the  combined  effects  of 
pressure  difference  and  wall  loss  on  the  local  atom  concentrations  between  two  PMT  positions  along 
the  flow  tube. 

For  N-atom  LIF  calibration,  flow  tube  measurements  were  made  without  added  NO  flow.  For 
O-atom  LIF  calibration,  NO  was  added  just  below  the  titration  endpoint.  Table  1  lists  the  pressure 
ratios  for  the  different  PMT  locations  under  our  flow  conditions,  as  well  as  the  estimated  ratios  of 
atom  concentration  for  different  assumed  values  of  y  .  Pressure  variations  are  small,  amounting  to 
only  a  few  percent.  However,  wall  losses  can  potentially  change  the  atom  concentrations  by  several 
10’s  of  percent  along  the  flow  tube,  depending  on  the  loss  probabilities  assumed.  We  choose  the  value 
Yn  -  Yn  ~  1 x  1 0  '  in  our  estimates  of  wall  loss  and  show  later  that  this  produces  reasonable  agreement 

between  simulations  and  experiment.  Table  2  compares  the  measured  PMT  signal  ratios  in  flow  tube 
mode  to  the  pressure  and  wall  loss  corrected  PMT  signal  ratios. 

Table  1.  Calculated  pressure  and  concentration  ratios  in  flow  tube  configuration 

Pressure  N-atom  Concentration  Ratios  O-atom  Concentration  Ratios 


PMT  Ratl0S  1x10~5  5x10~5  lx10~*  yq  =  Ixl0~5  5x10~5  Ixl0~* 


1 

2 

1 

0.982 

1 

0.979 

1 

0.900 

1 

0.809 

1 

0.980 

1 

0.906 

1 

0.821 

3 

0.964 

0.959 

0.811 

0.658 

0.962 

0.822 

0.676 

4 

0.946 

0.940 

0.733 

0.537 

0.943 

0.747 

0.559 

Table  2.  Measured  and  corrected  PMT  signal  ratios  in  flow  tube  configuration 


N-atom  PMT  Signal  Ratios 

O-atom  PMT  Signal  Ratios 

PMT 

Measured 

Corrected 

Measured 

Corrected 

1 

2 

1 

1.06 

1 

1.10 

1 

1.21 

1 

1.26 

3 

2.90 

3.14 

3.53 

3.81 

4 

0.694 

0.780 

0.565 

0.633 
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Relative  atom  concentrations  during  diffusion  tube  experiments  can  now  be  assigned  by 
dividing  the  L1F  data  by  the  corrected  PMT  signal  ratios  of  Table  2.  Figures  4a  and  4b  show  these 
relative  atom  concentrations  as  a  function  of  added  NO  flow,  for  N  atoms  and  O  atoms,  respectively, 
where  the  data  have  been  adjusted  to  the  same  endpoint  of  0.0425  seem  added  NO  flow.  The 
corresponding  N-  and  O-atom  concentrations  at  the  side-arm  entrance  can  be  approximated  by  the 
expressions  cQ  =  2.1 1x10*' F^/0.0425  and  cN  =  2.1  lx  10~'  (0.0425  -  FWO)/0.0415 ,  where 
concentrations  are  in  mol  m'  and  NO  flows  are  in  seem. 


0.00  0.01  0.02  0.03  0.04  0.05  0.06 


ADDED  NO  FLOW,  seem 


Figure  4.  Relative  atom  concentrations  at  the  different  PMT  locations  during  diffusion  tube 
experiments:  a)  N  atoms,  and  b)  O  atoms.  The  error  bars  are  the  standard  deviations  of  the  collected 
PMT  signals. 
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Reactor  Modeling 

A  numerical  reaction-diffusion  model  was  developed  at  Princeton  to  try  to  reproduce  the 
curvatures  and  relative  magnitudes  of  the  experimental  N-  and  O-atom  concentrations  plotted  in  Fig. 
4,  and  the  corresponding  scaled  measurements  in  Fig.  3,  using  the  entrance  atom  concentrations 
deduced  above  and  various  surface  reaction  boundary  conditions.  The  side-arm  reactor  geometry  was 
modeled  as  tube  of  length  L  =  1 46.7  cm  and  uniform  radius  R  =  1 . 1  cm,  closed  at  one  end  by  a  disk 
normal  to  the  tube  axis.  The  governing  equation  describing  two-dimensional  cylindrical  diffusion- 
reaction  of  species  s  in  the  side-arm  tube  is 


yy  +  -  Y  (rjsr )  +  y  0« )  = 
at  r  dr  dz 


(2) 


where  ps  is  the  species  mass  density,  w,  represents  the  rate  of  species  production  due  to  gas-phase 
chemical  reactions,  and  jsr  and  jK  are  the  diffusive  mass  fluxes  in  the  radial  and  axial  directions 
relative  to  the  mass  average  velocity  of  the  mixture.  The  diffusive  mass  fluxes  in  the  radial  and  axial 
directions  are  related  to  the  corresponding  diffusion  velocities  by  jsr  =  p,V„  and  ja  =  psvs. . 

Concentration  profiles  were  obtained  by  the  simultaneous  solution  of  the  set  of  8  species 
continuity  equations  (Eq.  2)  for  N,  O,  N2,  O2,  NO,  N2O,  NO2,  and  O3.  Solution  of  the  momentum  and 
energy  equations  is  not  required  to  obtain  concentration  profiles  under  the  isothermal  and  isobaric 
conditions  of  the  dead-end  side-arm  reactor.  The  ideal  gas  law  was  used  as  the  equation  of  state.  The 
species  mass  flux  vector  is  given  by 


js  =  -{cM sDs  )V.r(  +^£(cMrD,)V*r 

P  r 


(3) 


with  M s  the  species  molar  mass,  xs  the  species  mole  fraction,  and  Ds  the  effective  binary  diffusion 
coefficient  for  species  s  in  the  gas  mixture;  c  and  p  are  the  total  molar  and  mass  densities  of  the  gas, 
respectively. 

Effective  diffusion  coefficients  in  the  multi-component  gas  mixtures  of  the  experiment  are 
computed  usin^  the  friction-weighted  self-consistent  effective  binary  diffusion  (SCEBD) 
approximation  2'33: 


D.  = 


1  - 


<0  j 


Y  r 

I  — 
l  £?*>„, 
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Different  choices  are  possible  for  the  weighting  factor  0)K  32-33 :  we  choose  (Os  =  xs  so  that  Eq.  4 
reduces  to  the  exact  binary  value  for  a  two-component  mixture.  The  binary  diffusion  coefficients  Dr  v 
for  all  species  pairs  are  computed  from 


3_  l24kBTyNAn 

16  PaflMy  Mrs 


(5) 


where  M r  (  =  M rM  J(M r  +  M s)  and  are  collision  integrals  compiled  by  Wright  et  al.34  35.  The 

collision  integrals  were  implemented  using  the  curve  fit  formulas  of  Gupta  et  al.'6,  in  a  manner 
identical  to  that  employed  in  the  widely  used  NASA  CFD  code,  DPLR.  Collision  integrals  for  ozone 
were  unavailable,  and  were  approximated  by  those  given  for  carbon  dioxide.  Table  3  lists  the  room- 
temperature  rate  coefficients  for  nine  gas  phase  reactions  included  in  our  model  involving  dissociated 
oxygen  and  nitrogen  species.  The  reverse  reaction  rates  are  negligibly  slow  under  our  experimental 
conditions  and  were  set  to  zero. 


Table  3.  Gas-phase  reactions 


Reaction 

Rate3 

References 

R1 

0  +  0  +  M  — ^  O2  +  M 

kt  =(3.1±1.6)xlO  33 

38 

R2 

O  +  O2  +  M  — >  O3  +  M 

k2  =(6.0±0.6)xl0  34 

20 

R3 

0  +  O3  — >  O2  +  O2 

k3  =(8.0±1.2)xl0  15 

20 

R4 

O  +  N+  M  — >  NO  +  M 

k4=(  1.0±0.5)xl(T32 

21 

R5 

O  +  NO  +  M  — >  N02  +  M 

k5  =(9.0±1.8)xl0'32 

20 

R6 

O  +  NO2  — >  NO  +  O2 

k6={  1.04±0.15)xl0  11 

20 

R7 

n  +  n  +  m->n2+m 

k7  =(4.4±2.2)xl0'33 

21 

R8 

n  +  no->n2+o 

ks  =(3.0±0.9)xl0'u 

20 

R9 

N  +  NO2  ->  N20  +  0 

k9  =  (l.2±0.6)xl0" 

20 

a  Bimolecular  rate  constants  are  in  units  of  cm'  molecule'1  s'1.  Termolecular  rate  constants  are  in  units 
of  cm6  molecule'2  s'1,  and  the  collision  partner  M  is  taken  as  Nj. 

The  production  of  species  by  surface-catalyzed  reactions  must  be  balanced  by  their  diffusive 
flux  at  the  wall: 


= 


(6) 
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For  reactant  species  the  production  rates  at  the  surface  can  be  expressed  as  a  loss  probability,  ys , 
times  the  surface  impingement  fluxes;  for  O  and  N  atoms  this  gives  the  production  rates  (negative  for 
loss) 


(7a) 


(7a) 


The  loss  probability  is  defined  as  the  fraction  of  impinging  reactant  flux  removed  permanently 
from  the  gas  phase.  It  is  not  a  fundamental  chemical  quantity,  rather,  it  reflects  the  total  efficiency  of 
all  operating  surface  reaction  pathways  that  remove  species  5  on  a  particular  surface  under  a  particular 
combination  of  temperature,  pressure,  and  gas  composition.  If  the  reactant  species  is  consumed  to 
produce  more  than  one  product  species,  branching  fractions,  fsr ,  can  be  defined  as  the  fraction  of 
impinging  reactant  removed  permanently  from  the  gas  phase  that  participates  informing  product 

species  r.  We  avoid  the  common  term  “recombination  coefficient”  for  ■'  because  this  term  implies  a 
reaction  product  that  is  the  parent  molecule  of  the  reactant  (i.e.,  O2  for  O-atoms),  which  is  clearly  not 
always  the  case  in  multi-reactant  mixtures. 

The  absolute  lower  and  upper  bounds  on  both  ys  and  ft  r  are  0  and  1 ,  but  this  range  may  be 

reduced  by  other  factors  such  as  the  availability  of  partner  reactants.  For  a  partially  dissociated 
mixture  of  oxygen  and  nitrogen  interacting  with  a  catalytic  surface  via  three  net  steady-state 
heterogeneous  reaction  pathways:  O  +  O  — >  O2,  N  +  N  — >  N2,  and  O  +  N  — >  NO,  element  conservation 
dictates  the  constraint: 


(8) 


In  the  simulations  that  follow,  we  demonstrate  the  effect  of  NO  surface  production  on  the 
predicted  O-  and  N-atom  concentration  profiles  along  the  diffusion  side-arm  tube  for  different  fixed 
N2/N/O  compositions  at  the  side-arm  opening.  To  achieve  this  we  choose  values  of  y0  and  yN  ,  and 

then  vary  nitric  oxide  surface  production.  When  y0p0wv0  /M 0  <  yN pN  JvN / M  N  ,  the  level  of  NO 
production  can  be  set  by  choosing  0  <  fo  a  <  1  together  with 


(9a) 
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When  YoPo.*vo/M o  >  YnPn.»vn  /M  N  ,  NO  production  is  determined  by  choosing  0  <  fN  N_  <  1  and 


fo.o2  —  ^  0  /n.N2  ){YhPn,wVN^ o/YoPo,kvo^ n  )  • 


(9b) 


Nitric  oxide  surface  production  is  respectively  eliminated  or  maximized  when  /0  0  =1  or  0  in  Eq. 
(9a),  or  when  fN  Ni  =1  or  0  in  Eq.  (9b). 

The  individual  production  rates  of  the  surface  reaction  products  are 
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Equation  (6),  together  with  expressions  (7a,b)  and  (lOa-c),  provides  the  species  boundary 
conditions  on  all  the  side-arm  surfaces  at  r  =  R  for  z  =  0  to  L  and  at  z  =  L  for  r  =  0  to  R  .  A  zero- 
flux  boundary  condition,  jsr  =0,  was  applied  along  the  tube  centerline  ( r  =  0  for  z  =  0  to  L).  The 
gas  composition  at  the  side-arm  entrance  was  specified  according  to  the  experimental  conditions  and 
taken  as  uniform  for  z  -  0  and  r  =  0  to  R . 

Equation  2  was  discretized  using  a  fourth-order  accurate  finite  difference  scheme  for  the 
spatial  derivatives  and  iterated  to  steady-state  using  a  third-order  accurate  Runge-Kutta  method.  Grid 
convergence  was  obtained  for  grids  with  60  x  10  grid  points  in  the  axial  and  radial  directions, 
respectively.  Steady-state  solutions  were  achieved  using  a  CFL  (Courant,  Freidricks,  Levy)  number  of 
1 .  The  CFL  number  is  defined  as  the  ratio  of  the  simulation  time  step,  At ,  to  the  minimum  physical 
diffusion  time  scale  in  the  flow,  A tdjff,  where  A tdlff  =  p{Az 2  +A r2)/ p  and  Azand  Ar  are  the  grid 

spacings  in  the  axial  and  radial  directions,  respectively.  The  simulations  that  follow  were  all 
computed  using  a  60  x  10  grid  and  a  CFL  number  of  1 . 


Simulation  Results 

Simulations  were  computed  for  T  =  298  K,  P  =  60  Pa,  and  various  O-atom  and  N-atom 
concentrations  at  the  side-arm  entrance  corresponding  to  NO  titration  flows  in  the  range  0.0  to  0.0425 
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seem,  with  the  balance  assumed  to  be  molecular  nitrogen.  The  concentrations  of  all  other  species  were 
set  to  zero  at  the  entrance. 


ADDED  NO  FLOW,  seem 


Figure  5.  Model  simulations  for  y0  =  yN  =  lx  10  :  a)  relative  atom  concentration  profiles  along  the 

diffusion  tube  side-arm  at  an  added  NO  flow  rate  of  0.0215  seem,  with  or  without  surface  NO 
production;  b)  relative  atom  concentrations  at  different  PMT  locations  as  a  function  of  the  flow  of  added 
NO,  with  maximum  NO  surface  production. 


Figure  5a  shows  relative  atom  concentration  profiles  along  the  dead  end  diffusion  tube 
computed  for  an  NO  titration  flow  of  0.02 1 5  seem  and  the  loss  probabilities  yQ  =  yN  =  1  x  1 0'5 .  The 

solid  lines  show  the  decay  of  N-  and  O-atom  concentrations  along  the  diffusion  tube  when  no  surface 
NO  production  is  allowed.  Though  the  loss  probabilities  are  set  equal  for  N  and  O  atoms  in  this 
simulation,  N  atoms  are  lost  more  rapidly  down  the  side-arm  tube  than  O  atoms  because  the  lighter  N 
atoms  collide  with  the  surface  more  often  and  have  a  lower  diffusion  coefficient  in  molecular  nitrogen 
(about  522  vs.  628  cm2  s’1).  The  corresponding  profiles  when  NO  surface  production  is  maximized 
(dotted  lines),  show  that  the  loss  rate  of  N  atoms  increases  and  the  loss  rate  of  O  atoms  decreases. 

With  surface  NO  production  eliminated,  the  same  concentration  profiles  (Fig.  5a,  solid  lines)  are 
produced  whether  or  not  gas-phase  chemistry  is  included  in  the  computations.  With  surface  NO 
production  maximized  and  the  N  +  NO  gas  phase  reaction  eliminated,  computed  O-  and  N-atom 
profiles  again  overlay  the  profiles  for  no  NO  surface  production  (solid  lines).  Taken  together,  these 
modelling  results  support  the  mechanism  proposed  by  Copeland  et  al.11,  that  requires  both  surface 
production  of  NO  and  the  N  +  NO  gas-phase  reaction  to  explain  the  experimental  measurements. 
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In  Fig.  5b,  we  show  the  relative  N-atom  (solid  lines)  and  O-atom  (dotted  lines)  concentrations 
computed  at  the  different  PMT  locations  as  a  function  of  added  NO  flow.  These  computations  were 
made  assuming  maximum  NO  surface  production  and  yQ  =  yN  =  1  xlO'5 .  Comparison  with  the 

experimental  data  of  Fig.  4  shows  generally  good  agreement,  keeping  in  mind  that  the  initial  rise  of  N- 
atoms  with  NO  addition  is  not  included  in  the  reaction-diffusion  model.  Relative  atom  concentrations 
under  un-titrated  and  fully-titrated  conditions  are  independent  of  surface  or  gas-phase  NO  chemistry. 
Model  predictions  using  y0  =  yN  =  1  xlO  s  reproduce  measured  relative  atom  concentrations  at  0.0 
and  0.0425  seem  added  NO  fairly  well.  Much  poorer  agreement  is  obtained  for  higher  loss 
probabilities  like  yQ  =  yN  =5x10  '  which  predict  atom  losses  much  greater  than  those  observed 
experimentally. 


N-atom 


O-atom 


Figure  6.  Experimental  (symbols)  and  computed  (lines)  scaled  atom  concentrations  at  the  four  PMT 
locations  as  functions  of  added  NO  flow.  Computations  with  yQ  =  yN  =  1  XlO  and  a),  b)  maximum,  c), 

d)  intermediate,  and  e),  f)  minor  NO  surface  production.  Dotted  lines  show  a  linear  dependence  for 
reference. 


Figures  6a-f  compare  experimental  scaled  N-atom  and  O-atom  concentrations  to  numerical 
predictions  assuming  y0  =  yN  -  1x10  '5  and  different  levels  of  NO  surface  production.  The 
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experimental  data  have  been  scaled  to  align  (approximately)  the  experimental  and  computational 
endpoints  at  0.0  seem  (no  O  atoms)  and  0.0425  seem  (no  N  atoms). 

Figures  6a  and  6b  show  that  the  model  computations  with  maximum  NO  surface  production 
reproduce  the  concave  and  convex  curvatures  observed  in  the  experimental  diffusion  tube  data  for 
both  N  and  O  atoms.  However,  the  predicted  differences  in  curvature  between  four  PMT  locations  are 
much  larger  than  observed,  suggesting  that  the  maximum  surface  NO  production  condition  is  too 
extreme. 

A  similar  comparison  is  shown  in  Figs.  6c  and  6d,  for  calculations  that  enforce  a  more 
moderate  NO  surface  production  condition,  by  setting  either  fGO  =  0.5  in  Eq.  9a  or  fNN  =  0.5  in 

Eq.  9b.  The  computed  profiles  now  show  a  tighter  spread  in  curvature,  more  in  line  with  experimental 
observation,  especially  for  oxygen  atoms.  While  the  quality  of  agreement  depends  somewhat  on  the 
scaling  used  to  align  the  experimental  and  computed  endpoints,  the  general  observations  of  the 
experiment  are  well  reproduced. 

Figures  6e  and  6f  show  computational  results  when  NO  surface  production  is  a  minor  pathway 
(either  fQj0  =  0.9  in  Eq.  9a  or  fN  Ni  -  0.9  in  Eq.  9b.)  The  profiles  show  much  less  curvature, 

approaching  the  linear  dependence  indicted  by  the  dotted  line,  and  compress  the  differences  in 
curvature  between  the  four  PMT  locations  further.  These  simulations  do  not  reproduce  the  experiment 
as  well  as  those  of  Figs  6c  and  6d,  indicating  that  NO  surface  production  is  not  a  minor  pathway  in  our 
experiment,  but  rather  comparable  in  magnitude  to  the  O  +  O  and  N  +  N  surface  recombination 
pathways. 

The  reaction-diffusion  model  used  to  simulate  the  diffusion  tube  experiments  is  only 
constrained  by  element  conservation  and  does  not  include  finite  rate  surface  chemistry.  This  type  of 
model  is  often  used  for  parametric  studies  seeking  to  bound  catalytic  effects^9,  but  lacks  physics-based 
kinetic  formulations  that  describe  intermediate  catalytic  behavior  under  transient  thermal  and  flow 
environments.  A  more  attractive  approach  for  including  the  NO  surface  formation  in  CFD 
computations  is  with  a  finite-rate  surface  chemistry  model  as  implemented  by  Kurotaki  and  others^' 
6,8,9 ;  i.e.,  a  model  that  incorporates  kinetic  mechanisms  like  adsorption,  thermal  desorption,  Eley- 
Rideal  recombination  and  Langmuir-Hinschelwood  recombination.  Unfortunately,  these  models 
contain  large  numbers  of  numerical  parameters  that  must  be  chosen  by  theory  and/or  adjusted  to 
reproduce  experiment  data.  More  experimental  data  are  required,  at  higher  temperatures  and  over  a 
larger  range  of  gas  pressures  and  compositions,  to  enable  the  construction  of  a  suitable  finite-rate 
chemistry  model. 

Nevertheless,  the  experimental  results  presented  here  support  the  surface-catalyzed  production 
of  nitric  oxide  from  the  recombination  of  atomic  oxygen  and  atomic  nitrogen  on  quartz  surfaces. 
Computational  simulations  of  the  side-arm  reactor  with  a  multi-species  reaction-diffusion  model, 
incorporating  NO  surface  production  and  gas  phase  reaction  of  NO  and  N  atoms,  qualitatively 
reproduce  the  experimental  observations,  and  suggest  that  NO  surface  formation  is  of  comparable 
magnitude  to  the  O  +  O  and  N  +  N  surface  recombination  reactions. 
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SURFACE  PRODUCTION  OF  NITRIC  OXIDE  AT  ELEVATED  TEMPERATURES: 
EXPERIMENTS  USING  DIRECT  DETECTION  OF  NO 


The  experimental  approach  we  used  to  infer  the  production  of  NO  on  silica  surfaces  at  room 
temperature  can  in  principle  be  extended  to  elevated  temperatures.  However,  at  elevated  temperatures 
such  experiments  encounter  significant  complications.  The  primary  difficulty  is  the  rapidly  increasing 
intensity  of  the  blackbody  radiation  with  increasing  furnace  temperature,  which  creates  a  massive 
signal  background  in  the  745  nm  and  845  nm  observation  windows  used  to  detect  N  and  O  atoms  by 
2-photon  LIF.  The  atom  number  densities  in  the  diffusion  tube  at  high  temperatures  are  lower  than  at 
room  temperature  because  of  the  ideal  gas  law,  NHigh  =  N Ulw[TUw /THigh  ]  at  constant  pressure.  In 

addition,  the  heterogeneous  loss  of  atoms  on  hot  silica  walls  generally  increases  with  temperature  and 
the  LIF  scheme  for  O-atoms  probes  the  population  of  the  ground  state  with  the  J  =  2  fine  structure 
component,  which  also  decreases  with  rising  temperature/'  1  Thus,  extracting  small  LIF  signals  from 
the  overwhelming  thermal  radiation  background  at  high  temperature  becomes  extremely  challenging. 

We  have  developed  a  new  experimental  method  for  the  investigations  of  NO  surface 
production  at  elevated  temperatures,  based  on  the  direct  detection  of  NO  by  single-photon  LIF 
combined  with  photon  counting.  Surface-produced  NO  molecules  are  rapidly  lost  in  the  gas-phase  by 
the  reaction  NO  +  N  — »  N2  +  O;  however,  the  steady-state  NO  densities  resulting  from  the  balance 
between  the  surface  production  and  phase-phase  reaction  loss  are  still  well  within  the  reach  of 
sensitive  single-photon  LIF  detection.  The  LIF  from  NO  is  observed  in  the  ultraviolet  spectral  region, 
thus  thermal  radiation  background  from  the  hot  furnace  has  a  much  smaller  effect  than  in  the  near- 
infrared.  As  an  added  advantage,  NO  signal  detection  actually  becomes  easier  at  elevated 
temperatures,  because  the  NO  surface  production  is  expected  to  become  more  efficient,  while  the  rate 
coefficient  for  NO  removal  by  NO  +  N  — »  N2  +  O  decreases  with  increasing  temperature  ". 

The  experimental  setup  is  very  similar  to  the  one  described  previously  and  shown  in  Fig.  1 , 
with  the  following  differences.  For  NO  detection,  we  use  two  Hamamatsu  R928  PMTs  mounted  at 
different  distances  along  the  diffusion  tube.  Compared  to  the  R636  PMTs  used  for  N-  and  O-atom  LIF 
detection,  the  R928  PMTs  have  a  higher  sensitivity  in  the  ultraviolet  and  a  greater  internal  gain,  which 
makes  them  more  suitable  for  photon-counting  detection  of  the  248  nm  LIF  from  NO.  Instead  of  the 
gated  integrators,  we  use  two  SRS  SR400  photon  counters  to  record  the  PMT  signals,  which  allow  us 
to  work  with  fluorescence  signals  lower  than  1  photon  per  laser  pulse.  LIF  was  induced  by  laser 
output  with  wavelengths  near  226  nm,  corresponding  to  the  A—X  (0,  0)  band  of  NO,  and  fluorescence 
was  detected  near  248  nm  (A-X  (0,  1)  band)  using  PMTs  filtered  by  248  nm  interference  filters  with  a 
FWHM  of  10  nm.  To  avoid  saturation  effects  that  would  cause  nonlinearities  of  the  measured  LIF 
signal  with  respect  to  the  excitation  energy,  the  excitation  laser  beam  was  attenuated  using  a  Schott 
colored  glass  filter  (2  mm-thick  UG-5). 

The  main  difficulty  we  faced  using  this  new  approach  was  to  discriminate  between  the  desired 
signal  -  that  of  the  surface-produced  NO  -  and  background  signals  from  NO  present  in  the  system  from 
other  sources.  We  have  addressed  this  problem  by  employing  isotopically-enriched  titration  gases  as 
described  below. 

The  first  complication  arises  from  the  fact  that  the  titration  agent  is  the  same  species  as  the  one 
whose  surface  production  is  investigated.  The  titration  endpoint  is  never  particularly  sharp,  in  part  due 
to  instabilities  in  the  gas  flows,  implying  that,  below  but  near  the  titration  endpoint,  the  surface- 
produced  NO  signal  will  be  superposed  with  the  emerging  signal  from  the  NO  titration  gas  that  was 
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not  completely  consumed  in  the  gas-phase  reaction  with  N  atoms.  This  makes  it  impossible  to  discern 
the  behavior  of  the  surface-produced  NO  signal  near  the  endpoint. 

To  bypass  this  problem,  we  used  NO  gas  labeled  by  the  15N  isotope,  as  the  titration  agent.  The 
main  flow  of  nitrogen  gas  passing  through  the  discharge  produces  14N  atoms  and  the  nN  atoms  from 
the  ,sNO  titration  gas  are  incorporated  into  l5Nl4N  molecules  during  the  titration  reaction.  Any 
residual  titration  gas  in  the  system  would  be  |SNO  and  any  surface-produced  nitric  oxide  molecules 
would  be  i4NO.  The  measured  rotationally-resolved  LIF  excitation  spectra  of  l5NO  and  l4NO  are 
different  and  it  is  easy  to  trace  the  evolution  of  LIF  signals  of  both  species  with  the  changes  in  the 
flow  rate  of  the  titration  gas.  (The  natural  ratio  4N/  5N  is  about  270  so  strictly  speaking  some  |SN  is 
generated  in  discharge  and  some  '^NO  could  be  produced  on  the  surface,  but  these  levels  would  be 
more  than  two  orders  of  magnitude  smaller  than  surface  produced  l4NO.) 

Gas-phase  exchange  reactions  between  isotopically-labelled  nitrogen  atoms  and  nitrogen 
molecule  or  nitric  oxide  molecules  would  complicate  the  interpretation  of  these  experiments. 
Fortunately  the  rates  for  these  exchange  reactions  turn  out  to  be  extremely  slow.  Back  and  Mui41'  have 
set  an  upper  limit  for  the  rate  constant  of  the  reaction 14  N+'sNl'N^l'N+l4Nl'iN  of  6.6xl0~16  cm* 
molecule  1  s  'at  a  temperature  of  1000°C.  The  exchange  reaction  rate  decreases  with  decreasing 
temperature  and  so  the  rate  constant  should  be  even  smaller  near  room  temperature.41 4"  Back  and 
Mui4''  also  set  the  upper  limit  for  the  exchange  reaction  14  N+|SNO— >|SN+l4NO  to  be  at  least  250 
times  slower  than  the  rate  of  the  14  N+15NO  — »  0+l4N|sN  reaction,  based  on  experiments  at  both 
room  temperature  and  380°C. 

A  second,  somewhat  perplexing  complication  arose  because  of  a  background  NO  LIF  signal 
observed  at  elevated  temperatures,  even  without  any  titration  gas  added  into  the  system,  indicating  that 
a  source  of  oxygen  atoms  other  than  the  NO  titration  agent  was  present  in  the  system.  This 
background  signal  was  only  observed  with  the  microwave  discharge  running  and  was  found  to 
increase  as  the  furnace  temperature  rose.  The  immediate  suspects  were  an  O2  impurity  in  the  N2  gas 
or  an  air  leak  into  the  system,  both  of  which  could  lead  to  the  production  of  O  atoms  in  the  discharge. 
We  replaced  the  standard  99.999%  purity  N2  gas  (<1  ppm  specified  O2  impurity)  by  a  special  N2  grade 
with  specified  O2  impurity  less  than  0.01  ppm,  however  the  background  NO  signal  remained  as  large 
as  before,  demonstrating  that  O2  impurity  in  the  N2  gas  was  not  the  source.  We  were  unable  to  find 
any  leaks  into  the  system  using  a  He  leak  detector,  however  the  base  pressure  of  the  evacuated  system 
once  isolated  from  the  pump  did  slowly  rise  either  from  out-gassing  of  the  system  walls  or  from  some 
undiscovered  leak. 

In  a  separate  experiment,  small  amounts  of  air  were  added  at  the  titration  inlet  at  flow  rates 
(~0.2  seem)  that  produced  pressure  rises  similar  to  the  natural  pressure  rises  in  the  isolated  evacuated 
system.  These  air  additions  increased  the  NO  background  signal  by  factors  2  to  4.  This  experiment 
shows  that  the  oxygen  impurity  need  not  be  introduced  upstream  of  the  discharge  to  act  as  a  source  of 
O-atoms.  It  is  possible,  though  unlikely,  that  some  O2  can  diffuse  against  the  bulk  flow  to  dissociate 
in  the  discharge.  Another  possibility  is  that  O  atoms  are  created  by  gas-phase  reactions  of  O2  with 
electronically  excited  nitrogen  molecules,  particularly  N2  (A)  and  N2  (B),  created  by  the  discharge  or 
by  N  +  N  gas-phase  recombination  in  the  post-discharge  region.4  44  , 

Being  unable  to  eliminate  the  background  NO  signal,  we  resorted  again  to  the  isotopic  labeling 
approach,  this  time  also  labeling  the  O  atoms  of  the  titration  NO  gas  by  the  ihO  isotope.  The  titration 
reaction  now  becomes 
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I80  +  i5NI4N. 


I4n+  i5Ni80  -*■ 


The  surface  reactions  involving  N  atoms  generated  in  the  discharge  and  O  atoms  originating  from  the 
titration  gas  now  become 


l80  +  (,4N-S)  —  l4Nl80  +  (S) 
i4N  +  (i80-S)  -*■  l4Nl80  +  (S). 


Therefore,  we  assured  that  the  detected  l  4NIH0  nitric  oxide  has  an  O  atom  that  is  provided  by  the 
titration  gas,  and  an  N  atom  that  is  not  provided  by  the  titration  gas.  The  atomic  purities  of  the 
isotopically  labeled  gases  we  used  are  (99%  N,  98%  O)  for  |SN180  and  98%+  N  for  l5Nl60.  The 
natural  ,60/,80  ratio  is  about  500. 

The  measurement  protocol  was  as  follows.  Upon  heating  the  furnace  to  a  target  temperature, 
we  first  measured  the  LIF  excitation  spectrum  with  discharge  running  at  full  power  and  without  added 
NO  titration  agent.  The  spectrum  was  collected  by  slowly  scanning  the  dye  laser  wavelength  while 
counting  photons  in  both  PMT  channels  for  a  number  (10-50)  of  laser  pulses.  The  counting  period 
was  selected  such  that  it  did  not  diminish  the  wavelength  resolution  of  the  spectrum.  For  all  elevated 
temperatures,  the  spectrum  showed  features  easily  assigned  to  14NI60,  which  constitute  the 
background  signal  discussed  in  the  text  above. 

Isotopically  labeled  NO  titration  gas  was  then  introduced  and  the  spectrum  was  re-measured  at 
different  flow  rates.  The  titration  flow  rate  was  increased  in  small  increments  until  the  endpoint  was 
passed  as  indicated  by  the  disappearance  of  the  surface-produced  NO  signal  and  the  appearance  of  the 
LIF  signal  from  the  NO  species  introduced  by  titration. 

Finally,  we  performed  the  calibration  step:  the  discharge  was  turned  off,  NO  was  introduced  in 
the  system  with  a  known  flown  rate,  and  its  spectrum  was  measured.  The  calibration  has  two  roles. 
First,  it  serves  to  cross-calibrate  the  two  PMTs,  which  have  different  sensitivity  factors.  Those  factors 
can  be  determined  by  comparing  the  PMT  signals  in  the  calibration  run,  when  both  PMTs  detect  LIF 
corresponding  to  the  same  NO  number  density.  Second,  the  LIF  signal  strength  from  a  known  NO 
number  density  can  be  used  to  extract  the  absolute  number  density  of  the  surface-produced  NO  by 
comparison. 

Ideally  l4Nl80  would  be  used  for  this  calibration  step,  since  this  is  the  species  produced  by 
surface  reaction  in  the  experiments  with  the  1  N|H0  titration  gas.  Unfortunately  l4N180  was  not 
commercially  available  during  the  timeframe  of  this  study.  Therefore  we  use  experiments  with  |SNI80 
titration  gas  to  determine  the  general  trends  of  the  (unambiguously)  surface-produced  l4N180  number 
density  as  a  function  of  titration  gas  flow  rate.  Absolute  values  of  the  surface-produced  NO  number 
density  at  different  elevated  temperatures  were  determined  separately  from  '^Nl60  titration 
experiments  that  generated  i4N'  O.  The  LIF  signals  collected  during  these  experiments  were  corrected 
by  subtracting  the  background  l4Nl60  signals  discussed  above  and  then  comparing  the  signal  strength 
to  calibration  measurements  using  standard  l4Nl60  gas  (10%  NO  balance  N2).  This  procedure  adds 
some  uncertainty  to  the  results. 
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Figure  7.  LIF  spectra  at  1000  K  obtained  using  l5Nl60  titration.  The  spectra  were  measured  at  a 
distance  of  100  cm  down  the  diffusion  tube  and  represent  the  sum  photon  count  for  40  laser  pulses.  Top: 
spectrum  measured  below  the  titration  endpoint,  showing  lines  of  l4Nl6().  The  lines  that  were  used  to 
calculate  the  changes  of  LIF  intensities  with  the  titration  flow  rate  (Fig.  8)  are  labeled  by  arrows. 
Bottom:  spectrum  measured  above  the  titration  endpoint,  showing  lines  of  15N160. 


To  minimize  the  data  acquisition  time,  we  scanned  the  dye  laser  wavelength  over  relatively 
narrow  spectral  regions  that  contain  several  well-separated  lines  of  each  of  the  NO  species  of  interest. 
Figure  7  shows  spectra  collected  in  the  experiments  at  1000  K  using  15NI60  as  the  titration  agent. 
Below  the  endpoint  (observed  at  a  l<;Nl60  flow  rate  of  -0.57  seem),  the  spectrum  shows  practically 
only  i4Ni60  features;  the  very  weak  |SN160  features  in  the  226.30-226.32  nm  range  may  be  due  to 
residual  15NI60  left  over  in  the  diffusion  tube  from  the  previous  run  above  the  endpoint,  or  partly  due 
to  l5NlfD  naturally  present  in  the  14Nlf>0  cylinder.  The  spectrum  shows  raw  photon  count,  and  it 
illustrates  the  size  of  the  background  blackbody  radiation  signal  relative  to  the  LIF  signal.  Above  the 
endpoint,  only  strong  |fNl60  features  are  observed  (bottom  spectrum  at  Fig.  7).  For  either  species, 
there  are  several  lines  that  are  well  separated  from  the  lines  of  the  other  species  and  can  be  used 
unambiguously  for  species  identification. 

After  normalizing  the  spectrum  by  dividing  the  signal  with  the  measured  laser  pulse  energy,  4 
well-separated  l4Nl60  lines  were  fitted  with  Gaussians  and  the  sum  area  of  the  4  Gaussians  was  used 
to  trace  the  dependence  of  the  14NI60  LIF  signal  on  the  titration  gas  flow  rate.  This  dependence  is 
shown  in  Figure  8,  in  which  absolutely  calibrated  signals  for  both  PMTs  are  plotted.  Starting  from  a 
non-zero  value  at  zero  lsNlf,0  flow,  the  HNI60  number  density  increases  up  to  about  0.45  seem,  and 
then  rapidly  falls  towards  zero  as  the  flow  rate  approaches  the  endpoint  value. 
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Figure  8.  LIF  signal  intensity  for  14N160  as  a  function  of  the  15N160  titration  flow  rate.  Measurements 
were  performed  at  1200  K  for  two  different  distances  down  the  diffusion  tube  (as  indicated  in  the  legend). 
The  plotted  values  represent  integrated  photon  counts  for  several  selected  lines,  as  indicated  in  Fig.  7. 

The  signals  have  been  absolutely  calibrated  using  14N160  to  obtain  the  corresponding  number  densities 
(see  text). 


The  number  density  values  obtained  by  14NlhO  calibration  are  in  the  1011  cm  3  range,  which  is 
too  high  to  be  attributed  to  three-body  N  +  O  recombination  at  the  low  pressure  and  high  temperature 
used  in  the  experiment.  Therefore,  they  have  to  be  produced  by  surface-catalyzed  processes.  The 
background  NO  LIF  signal  is  sizable  at  this  high  temperature,  and  its  behavior  as  a  function  of  the 
titration  flow  rate  must  be  understood  separately  before  it  is  subtracted  from  the  total  signal.  This  was 

K  C  1  O 

made  possible  by  the  N  O  titration  experiments. 

1  C  1  o 

Figure  9  shows  spectra  measured  at  1200  K  using  N  O  titration.  A  different  spectral  region 
is  selected  than  in  the  case  of  1  Nl60  titration.  The  situation  is  also  more  complicated  in  this  case, 
because  both  l4N160  and  ,4NlxO  features  are  observed  in  the  spectrum  for  all  titration  flow  rates  below 
the  endpoint.  However,  l4NlhO  and  ,5NI80  spectra  that  we  measured  separately  allow  us  to  assign  the 
observed  14N160  transitions  unambiguously  and,  by  exclusion,  those  transitions  that  belong  to  neither 
i4N'\)  nor  1  NlxO  are  assigned  to  l4N|vO.  The  line  assignments,  as  well  as  the  lines  used  to  trace  the 
dependences  of  l4Nl60  and  1  *N|hO  number  densities  on  the  titration  flow  rate  are  marked  in  Figure  9. 
The  bottom  panel  in  Figure  9  shows  the  spectrum  measured  above  the  endpoint,  showing  only  the 
strong  lines  of  the  i5N*  O  titration  agent. 

Figure  10  shows  the  dependences  of  the  ^N'^O  and  I4NihO  number  densities  on  the  l5N180 
titration  flow  rate.  The  signals  from  the  two  PMTs  have  been  cross-calibrated,  as  described  above. 
The  uNlbO  background  signal  decreases  with  the  flow  rate  in  an  approximately  linear  fashion,  and  a 
linear  Fit  to  this  dependence  would  yield  a  non-zero  value  at  the  titration  endpoint.  Near  the  endpoint 
the  l4N1*0  signal  abruptly  drops  to  zero.  The  decrease  in  the  background  l4Nl0O  signal  with 
increasing  titration  flow  is  not  surprising,  given  that  no  new  ,hO  atoms  are  added  in  the  system, 
whereas  the  number  density  of  l4N  atoms  decreases  with  the  addition  of  |SNI80  titration  flow.  The 
strong  deviation  from  linear  dependence  near  the  endpoint  and  the  abrupt  drop  to  zero  are  harder  to 
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Figure  9.  LIF  spectra  at  1200  K  obtained  using  ,5N,80  titration.  The  spectra  were  measured  at  a 
distance  of  100  cm  down  the  diffusion  tube  and  represent  the  sum  photon  count  for  20  laser  pulses.  Top: 
spectrum  measured  below  the  titration  endpoint,  showing  lines  of  f4N,80  (labeled  by  ‘18’)  and  ,4N160 
(labeled  by  ‘16’).  The  lines  that  were  used  to  calculate  the  changes  of  LIF  intensities  with  the  titration 
flow  rate  (Fig.  10)  have  labels  enclosed  in  squares.  Bottom:  spectrum  measured  above  the  titration 
endpoint,  showing  only  15N18Q  lines. 


Figure  10.  LIF  signal  intensities  for  14N18G  and  l4N,60  as  functions  of  the  l5Nl80  titration  flow  rate. 
Measurements  were  performed  at  1200  K  for  two  different  distances  down  the  diffusion  tube  (as 
indicated  in  the  legend).  The  plotted  values  represent  integrated  photon  counts  for  several  selected  lines, 
as  indicated  in  Fig.  9.  The  signals  have  been  normalized  with  respect  to  the  excitation  pulse  energy  and 
cross-calibrated  with  respect  to  the  PMT  detection  sensitivities  (see  text). 


23 


understand  without  a  detailed  model.  However,  we  note  that  the  probed  steady-state  4NI60  number 
density  is  determined  by  a  competition  between  the  surface-catalyzed  production  and  collisional 
removal  by  N  atoms;  rates  for  both  processes  increase  with  the  N-atom  number  density,  and  in  such  a 
situation  strong  nonlinearities  near  zero  N-atom  number  density  (endpoint)  are  not  surprising. 

The  ,4N180  signal  plotted  in  Figure  10  is  zero  at  zero  titration  flow,  and  it  increases  with  flow 
rate  with  a  noticeable  curvature.  Similar  to  the  data  in  Figure  8,  the  curve  exhibits  a  rather  broad  peak 
near  2/3  of  the  endpoint  flow  rate,  and  then  rather  abruptly  decreases  to  zero.  The  14N160  data  in  Fig. 
8,  which  include  the  flow  rate-dependent  background  signal  contribution,  show  no  significant 
curvature  below  0.4  seem,  but  we  note  that  the  experimental  uncertainty  is  significant  enough  that 
such  curvature  cannot  be  excluded.  The  general  shape  of  the  l4NlxO  signal  dependences  plotted  in 
Figure  10  provides  key  guidance  for  our  surface  chemistry  model.  The  fact  that  the  14NlgO  signal  is 
purely  titration  induce,  i.e.,  it  disappears  at  both  zero  titration  flow  and  titration  endpoint, 
unambiguously  proves  that 1  lNl80  is  a  result  of  a  reaction  between  IH0  atoms  introduced  by  titration 
and  i4N  atoms  produced  by  N2  dissociation  in  the  discharge. 

For  the  purpose  of  initial  analysis,  we  assume  that  the  total  NO  number  density,  such  as  data 
plotted  in  Fig.  8,  is  a  simple  sum  of  the  background  and  titration-induced  number  densities  for  any 
titration  flow  rate.  We  can  then  extract  the  titration-induced  number  density  using  N  O  titration 
data  as  guidance.  At  zero  titration  flow,  all  of  the  NO  number  density  plotted  in  Fig.  8  is  a  result  of 
the  background  production.  At,  e.g.,  0.4  seem  titration  flow  rate,  Fig.  10  indicates  that  this  initial 
background  number  density  is  reduced  by  about  50%  for  a  distance  of  29  cm  down  the  tube.  Thus, 
from  the  total  NO  signal  at  0.4  seem  we  subtract  50%  of  the  initial  (zero-flow)  signal  to  obtain  the 
titration-induced  number  density.  Such  analysis  was  done  for  all  investigated  temperatures  and  the 
results  are  plotted  in  Fig.  1 1 . 
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Figure  11.  Number  density  of  surface-produced  14N,60  as  a  function  of  temperature  for  two  positions 
down  the  diffusion  tube.  The  data  was  obtained  using  1?N160  titration,  with  a  titration  flow  rate 
corresponding  to  70%  of  the  endpoint.  The  background  signal  was  subtracted  to  obtain  the  14N160 
density  due  solely  to  the  O  atoms  introduced  by  titration  (see  text). 
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The  room-temperature  NO  signals  are  very  small  and  the  error  bars  in  Fig.  1 1  hence  include 
zero;  the  method  is  not  well  suited  for  low  temperatures  and  it  is  thus  complementary  to  the  O-  and  N- 
atom  detection  approach  described  in  the  previous  section.  The  NO  density  exhibits  a  relatively  rapid 
increase  with  temperature,  except  for  the  possible  drop  between  600  and  800  K. 

Because  of  the  experimental  complexity  of  these  experiments,  sufficient  time  did  not  remain  to 
complete  finite-rate  surface  chemistry  modeling  of  these  results.  However,  the  same  model  used  in 
the  previous  set  of  experiments  shows  qualitative  agreement  with  the  measured  behavior  of  the  surface 
generated  NO. 


Table  4  Gas-phase  reactions  with  temperature  dependence 


Reaction 

Rate8 

References 

R1 

o  +  o  +  m->o2  +  m 

jfc,  =  1.57x10  34 7°exp(900/7) 

38 

R2 

O  +  02  +  M  — >  Oi  +  M 

k2  =  3.9xl0''297~2exp(-140/7) 

R2b 

O3  +  M  — >  O  +  O2  +  M 

k2b  =  4.13xl0~l(,7^’exp(-l  1430/7) 

38 

R3 

O  +  O3  ->02+02 

*3  =  8.6x  10~l27°exp(-2090/T) 

38 

R4 

O  +  N+  M  — »  NO  +  M 

k4  =  1.8xl0~3l7~°5exp(0/7) 

21 

R5 

O  +  NO  +  M  — »  N02  +  M 

k5  =  4.2  x  1  0_337°  exp(940/  7) 

21 

R6 

O  +  N02  — >  NO  +  02 

k6  =  1.7x10  n7°exp(-300/7) 

21 

R7 

n  +  n  +  m->n2+m 

k-,  =  8.3x10_347°  exp(500/7) 

21 

R8 

n  +  no->n2+o 

kg  =  2.7  x  10‘"7°exp(0/7) 

21 

R9 

N  +  N02  ->  N20  +  O 

k9  =  5.8x  10"i27()  exp(220/7) 

20 

8  Bimolecular  rate  constants  are  in  units  of  cm'  molecule  1  s'1.  Termolecular  rate  constants  are  in  units 
of  cm6  molecule  2  s'1,  and  the  collision  partner  M  is  taken  as  N2. 


The  temperature-dependent  gas-phase  reaction  rates  of  Table  4  were  used  for  these 
calculations.  The  general  form  for  these  reaction  rates  is  taken  as  =  AT”  exp (-(E / R)/T) ,  with 
specified  constants  A,  n  ,  and  El  R ,  which  can  sometimes  be  zero.  Because  ozone  thermally 
dissociates  at  appreciable  rates  above  -500  K,  one  more  equation  (R2b)  was  added  to  the  gas-phase 
reaction  set. 

In  the  simulations  shown  below,  the  entire  diffusion  tube  is  modeled  as  isothermal,  the  total 
pressure  is  58.5  Pa,  titration  endpoint  is  0.6  seem  in  a  73  seem  N2  flow,  and  the  branching  fractions  are 
specified  as  fN  Ni  =  /oc,  =  0.5. 
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Figure  12.  a)  Number  density  of  surface-produced  NO  as  a  function  of  distance  from  the  side-arm 
entrance  for  three  different  temperatures  and  two  different  reaction  efficiencies;  b)  number  density 
29  cm  and  100  cm  from  the  side-arm  entrance  as  a  function  of  temperature  for  two  different  reaction 
efficiencies.  The  dissociated  nitrogen  flow  was  2/3  titrated. 


Figure  13.  Number  density  of  surface-produced  NO  as  a  function  of  NO  added  at  the  titration  port,  at 
distances  of  29  cm  and  100  cm  from  the  side-arm  entrance. 


Figures  1 2a  and  1 2b  show  the  model  predictions  for  NO  number  density  as  a  function  to 
temperature,  distance  from  the  side-arm  entrance  and  reaction  efficiency.  The  NO  number  density 
rises  with  increasing  temperature  and  increasing  reaction  efficiency,  and  is  relatively  constant  along 
the  side-arm  tube  beyond  -20  cm  from  the  tube  entrance. 


26 


The  predicted  NO  number  densities  are  on  the  order  of  101 1  cm  3  for  a  reaction  efficiency  of 
10'4,  in  excellent  agreement  with  experiment,  considering  the  many  reaction  and  diffusion  parameters 
that  must  be  input  into  the  program.  We  note  that  computations  made  with  or  without  reaction  R4 
over  our  entire  temperature  range  produce  essentially  identical  numerical  results,  so  the  gas-phase 
production  of  NO  by  O  +  N  recombination  is  negligible  compared  to  the  surface  production  of  NO,  as 
stated  earlier. 

The  rise  of  NO  number  density  observed  experimentally  (Fig.  1 1 )  is  about  a  factor  of  four  from 
298  to  1200  K,  while  it  is  only  about  a  factor  of  two  over  this  temperature  range  in  model  predictions 
for  fixed  reaction  efficiency  (Fig.  12b).  However,  the  reaction  efficiency  for  atom  loss  on  quartz  is 
expected  to  increase  with  temperature,  so  the  agreement  between  computations  and  experiment  could 
easily  be  improved  by  adjusting  the  input  reaction  efficiency  as  a  function  of  temperature. 

Figure  13  shows  the  predicted  variation  of  the  surface-produced  NO,  29  cm  and  100  cm  from 
the  side-arm  entrance,  as  a  function  of  NO  added  at  the  titration  point.  The  calculations  are  for  a 
temperature  of  1200  K  assuming  a  reaction  efficiency  of  10  4.  Figure  13  can  be  compared  to  the 
experimental  '4N|hO  data  of  Fig.  10.  The  computations  do  a  good  job  of  reproducing  the  asymmetry 
of  the  NO  number  density  variation  with  titration  gas  addition,  rising  more  slowly  with  the 
introduction  of  NO  and  falling  more  rapidly  near  the  titration  end-point.  The  model  predicts  that  NO 
concentrations  at  100  cm  should  be  equal  to  or  greater  than  at  29  cm.  Experimentally  the  number 
density  is  usually  higher  at  29  cm  than  at  100  cm  (Fig.  10  and  1 1 ),  but  the  error  bars  on  the  measured 
number  densities  are  substantial  and  overlap  (Fig.  1 1). 

It  was  not  possible  implement  the  non-isothermal  tube  wall  temperature  into  the  computations 
at  this  time,  nor  to  do  a  sensitivity  study  on  the  gas-phase  reaction  coefficients,  in  particular  reaction 
R8.  However,  both  will  be  performed  and  the  conclusions  will  be  included  in  a  manuscript  currently 
under  preparation  describing  these  experiments. 


SUMMARY 

During  this  project  we  have  conclusively  demonstrated  the  importance  of  the  O  +  N  surface 
reaction  for  catalytic  atom  recombination  processes  on  silica  surfaces  in  dissociated  nitrogen-oxygen 
mixtures.  This  was  accomplished  using  a  diffusion  tube  side-arm  reactor,  a  microwave  discharge  to 
generate  atomic  nitrogen,  and  an  NO  titration  technique  to  generate  different  mixtures  of  O  and  N 
atoms,  together  with  two  distinct  species  diagnostic  approaches. 

In  the  first  approach,  two-photon  LIF  of  atomic  oxygen  and  atomic  nitrogen  were  used  to  show 
how  the  relative  atom  concentrations  varied  with  position  along  the  diffusion  tube  as  a  function  of  O/N 
ratio  at  the  diffusion  tube  inlet.  The  experimental  result  confirmed  and  extended  earlier  experiments 
by  Copeland  et  al."  and  were  consistent  with  a  reaction-diffusion  model  of  the  diffusion  tube  in  which 
O  +  N  surface  reaction  plays  a  comparable  role  to  O  +  O  and  N  +  N. 

In  the  second  approach,  single-photon  LIF  was  used  together  with  isotopically-labeled  titration 
gases  to  detect  surface-produced  NO  directly  at  two  different  positions  in  the  diffusion  tube,  as  a 
function  of  O/N  ratio  at  the  diffusion  tube  inlet.  These  experiments  independently  confirm  the 
production  of  NO  by  surface  reactions  on  quartz,  from  room  temperature  up  to  1 200  K.  The  same 
reaction-diffusion  model  used  to  interpret  the  previous  room-temperature  experiments  can  also 
reproduce  most  features  of  the  new  experimental  data  at  higher  temperatures. 
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